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Proxemy Research is under contract to NASA to perform science research on the effects of 
topography on geologic mass movements. The current project concentrates on the sensitivdy of 
lahars and debris flows to various scales of topographic variability. The following report constitutes 
delivery of Milestone Event #12 under NASA contract NAS5-99182. 

TITLE: 2000 Annual Report; ‘Topographic Effects on Geologic Mass Movements', PRI 

Tech. Rep. # TOPO-00-012. 

AUTHOR: Stephen M. Baloga 


1. Introduction 

Funding for this project began in April 1999. Since the project’s inception, we twve 
completed all milestones established in contract NAS5-99182. Monthly progress reports have beeti 
submitted, detailing work accomplished over each time period. This final report documents an 
summarizes the results of the entire project, including conclusions based on our findings and 
recommendations for further investigation. 

Proxemy’s efforts towards fulfilling these tasks has been divided into two areas, theoretical 
investigation and applications to actual mass movements. The theoretical investigation focused on 
determining differences in flow thicknesses and advance rates when the mass movements flow over 
topography with variability and roughness at different scales. Here we have used both steady-state 
and time-dependent models. We have also applied these models to actual lahars at Mt Ruapehu 
New Zealand. The applications have also included a major collaborative effort with Dr. Mark 
Bulmer (CEPS/NASM) to collect and analyze appropriate data from the recent debns flows m 

Madison County, Virginia. 

The research progress for each of the monthly milestones is summarized below. Details of 
the data products and scientific analyses appear as Appendices. In conclusion, we sumtnanze the 
scientific progress achieved by this effort and identify promising areas of future research. 

2. Summary of Milestone 1 

The initial efforts were focused on background investigations of well-documented lahars and 
similar mass movements at Mt. Ruapehu and Mount St. Helens. The goal was to begin compiling 
data for evaluating the sensitivity of mass movements to topographic variations and identify where 
the data sets were incomplete. 


3. 


Summary of Milestones 2, 3, and 4 
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These milestones addressed a theoretical estimation of the srasitivity of rapid m^ 
movements to variations in topography. A simple quantitative model was denved for 1^^ 
nroDagation assuming an empirical flow rate law, local and global conservation of flow, constan 
LffiLnt of resistant, and steady-state conditions. Such restrictions make the flow respond only 
locallv to topography variability. This elementary model was applied to a number of simuM 
tono Japhies^wSi identical average dimensions (elevation loss of 10 km over a 30 km length) bu 
JdL varying local slopes. The topography was simulated by 

varyii i^des and cycles to a uniformly sloping plane. The results show that both thicknesses 
front arrival tiLs for identical average topographies are highly dependent upon locd 
ZLlZZ lZits In particular, the presence of irregularities in the slope can sigmficantly slow 
the iS^ Flow front arrival times were seen to increase with increasing amplitudes and increasing 

considerably greater than that expected by increased travel distances along 
the curve of the roughened topography. These results show the importance of fine-scale topographic 
LTalrate prediction of hazards from lahars, debris flows, and similar mass movements. 

4. Siiminary of Milestones 5, 6, and 7 

A new time-dependent model for lahar transport was derived based on local global 
conservation of volume for a commonly used empirical flow rate. T^e f 

the thickness profile contain a term for cumulative topographic influence. T p y 

unrecoenized cumulative effect arises from three factors: nonlineanty m the local r^te, 

topogrSiic changes along the path of the flow, and time-dependence in the longitudinal thiclmess 
profit Vs model was compared to the steady-state model using identical underlying topographies. 

tpo^aW and the arrival times were shown to be dramatically different when tirn^ 
dependence was included. This, in turn, amplifies the need for topographic data at scales compar 
to the thickness of the lahar or debris flow. 

Proxemy Research personnel also began supporting an ongoing field 
1995 debris flows at Madison County, Virginia, by Dr. Mark Bulmer <CEPS/NAS^^ f 
this period were directed toward field measurements of topography using GPS ^d laser raging. 
ProxVy personnel were active in the acquisition of longitudinal and transverse channel profiles for 

one of the primary debris flows that occurred in 1995. 


5. Summary of Milestones 8, 9 and 10 

As part of a general field investigation of the Generals Landslide m Madison County, 
Virginia, a high-resolution topographic data set was and compiled, including f 
by CEfS and Proxemy Research. The initial results from this effort were presented at foe East Coas 
Volcanologists’ Meeting at CEPS during November, 1999. Results with 
appUcVons were written up in an abstract submitted to the Lunar and Planetary Science Conf e^e 
(Appendix 1 .A). Preliminary results from the Madison debris flow mvestiga ion P 
during a poster session at the Lunar and Planetary Science Conference (see Appendix 2). 
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6. Summary of Milestone Event 11 

Research activities focused on mass movements of the Cascade volcanoes. Topogr^hic data 
for Mount St. Helens and Mount Rainier were obtained ftom USGS DEMS^ A commercid 

CTopol”) was used to characteiire the slopes at Mount Rainier, wWe cross-sectioM ^eas 
^rlnber of historic mass movements were obtained from a USGS Professional Paper (#1547). 

Work also continued on the Madison County debris flow field site. The initial flow 
modeling results had proven unsatisfactory, as large, unexplained vanations m rheology appear^ 
to exist Numerous meetings were held to develop interpretations of the field data conipatible m th 
the context of a theoretical model. Agreement was reached on all 
alternatives for the dimensions of the flow. We now consider the compiled 

Landslide to be state^of-the art for a debris flow that was largely unobserved. This data set appears 
as Appendix 1 . The modeling efforts were reapplied with satisfactory results. A dra manuscnp 
for publication in Geophysical Research Letters was also begun (see Appendix 4). 


7. Summary of Milestone Event 12 

Significant effort was devoted toward the documentation of the data set and field 
interpretations for the Generals Landslide. Similar efforts were directed toward obtaining similar 
comprehensive data sets for mass movements at Mt. Ranier and Mt. St. Helens from the literature 
and other sources. An abstract for follow-on work was submitted to the Amencan Geophysica 
Union for the Spring 2000 meeting (see Appendix 3). Contributions to the draft GRL manuscnpt 
were made and effort was devoted to the final documentation of this project. 


8. Conclusions 

The final year of this project has been productive. Theoretical models have b^n developed 
for application to lahats, debris flows, and other geologic mass movements. These models have b^n 
formdated to be used in conjunction with field and remote sensing data. Using these tn»dds, we 
have investigated the effects of topography at a variety of scales on the emplacement 
flow profiles of lahars and debris flows. We have applied one of these models to a small wel - 
documented debris flow and have interpreted the results m the context of a general field 
investigation. 

We have shown that topographic resolution can significantly affect emplacement velocity, 
flow thicknesses, and estimates of flow rheology. We have also demonstrated that using mteipolated 
slope values derived from regression fits for high resolution data yields rnodehng results more 
reliable than those obtained with either low or high resolution topographic data. 

One future activity will be to apply the time-dependent model to the data sets for the 
Generals Landslide and similar mass movement at Mt Rainier and Mt St Helens. This effort should 
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lead to a gtoater understanding of the time-dependent effects of topography on flow emplacement 
and provide improved assessments of hazards. 

This research has produced two abstracts, 

—a. advances in undemtanding the sensitivity of 
rapid mass-movements to variations in topography. 
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APPENDIX 1; 

Investigation Methodology and Compiled Data for the Generals Debris Flow 

in Madison County, Virginia 


APPENDIX 1.A: 

LPSC Abstract 

APPENDIX 2: 
LPSC Poster 

APPENDIX 3: 
AGU Abstract 

APPENDIX 4: 


Draft GRL Manuscript 


APPENDIX 1* INVESTIGATION METHODOLOGY AND COMPILED DATA FOR 
THE GEN^S DEBRIS FLOW IN MADISON COUNTY, VIRGINIA 


Introduction: 

Detailed field studies of mass movements such as debris 
essential for interpreting data in the context of theoretical emplacment models. To tos end, 

Re“^h pe^Xihave participated in a GPS field survey of the Generals I^dsl.de tn Madtson 
Lunty Virginia, headed by Mark Bulmer of the Smithsonian Instttuhon s C^ter to “ 
Planed Studies (CEPS). The elevation of the main channel was survey^ in detml ‘E® 

“to debnl flow and information about the roughness of to " 
tramects were made across the chaimel at various stations along the flow path The ^ranse 
positions are overlayed on an air photo in Figure 1. The modeling data denved from 
ZSn^e. posLn, elevation, and slope; minimum and maximum width, area ‘h« 
is given in Table 1 . Additional information describing to transects is summanz 

Field survey: 

The main channel of the “General slide” was surveyed using a Trimble Total Station 4800 
which consists of two GPS receivers that work in tandem. This equipment uses c^er-ph 

at a known location during the survey. The location of this receiver was determined to accuracy 

of 1-2 m using a Trimble ProXRS GPS unit using satellite-based differential ^ . 

station was then used as the reference point for all subsequent surveys. The roving Tn 

conSfly conducting transects up the m^n 

m about lOOO m up from the lower end of to channel, until we encountered difficulties in satellite 

iCe due to scmening by to tree canopy. Beyond that point, the ™ 

cm. We successfrilly collected data up to the top of the mam channel, about 1600 m. The data point 

are shown in a scattergram (Figure 2). 

Slope algorithm: 

The elevation data obtained during the field survey were plotted as fimctions of downstream 
distance (Figure 3). A third-order polynomial regression was ton appli^to this data to denve a 
characterization function for the elevation, (Z) as a function of distance (X). 

Z = 5.617 X 10 * X^ + 2.544 x 10 '* X^ - 0.472968 X + 613.03 
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A gradient function was then obtained from the derivative of the elevation function: 
dZ/dX = 16.851 X 10 * + 3.088 x lO '* X - 0.472968 
Slopes were then calculated at each transect position from. 


0 = arctan (dZ/dx) 

Issues: Minimum and Maximum parameter values from field data 

Under actual field conditions, defining the shape of the active flow lobe from the remaining 
evidence proved a non-trivial task. Depth markers such as scars on trees and debns piles were used 
as indicators of both depth and position of the flow edge. However, in some cases the geomorphic 
evidence was absent or inconclusive. Air photos of the site taken immediately following the even 
proved useful supplements, but yielded direct information only on the channel and deposit, and not 
the active lobe itself. There was considerable evidence for superelevation, with markers indicating 
that flow depths on one bank may well have exceeded those on the opposite bank. Thus, it was not 
possible to absolutely constrain the dimensions and position of the active flow. Instead, minimum 
estimates (based on a horizontal flow surface at the lowest depth) and maximum estimates (assuming 
a “tilted” flow surface between depth indicators) were used to constrain the dimensions, with the 
actual values most likely falling between these two extremes. The transects, with minimum and 
maximum flow surfaces, are illustrated in Figure 4. For comparative purposes, the cross-sections 

are plotted at same scale in Figure 5. 


Cross-sectional area and thickness algorithm: 

To estimate the minimum cross-sectional area from the transect data (Figures 4 and 5), Each 
flow transect is divided into segments (Figure 6a). Each data point is the right edge of one se^ent 
and the left edge of the next (except for endpoints). The width of the segment is then the distance 
between the two data points that make the edges. The segments are treated as rectangles; average 
elevation of the bottom is calculated from an unweighted average of the two edge points. Depth is 
found by subtracting this elevation from the assumed elevation of the flow surface. The contnbution 
of each segment to area is then the segment width times this depth; these are summed to estimate 

cross-sectional area for the entire transect. 


The average thickness for the cross-section is found from the area divided by the width, 
effectively, a weighted mean. Thus the wider the segment, the more its depth is considered in the 
determination of the average. 

The method for determining maximum areas is basically the same, only the elevation of the 
top of each segment, instead of being the same horizontal flow level, is extrapolated from a linear 
equation, assuming a tilted straight-line superelevated flow surface. For each segment we calculate 
an average flow level between the two end points, assume the flow is “stepping up and thus flat 
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across the top, and find the thickness of the segment column from the difference between this 
elevation and the average elevation of the bottom-again dividing the flows into rectangles, but now 
each rectangle has its own top (Figure 6b). 


The Generic Empirical Rheology Model 

The initial model applied to the Generals debris flow data assumes steady-state flow volume 
conservation and a flow rate of the form: 


/ 

Q=uA = 

\ 


gsin0/i 


k 

A 


( 1 ) 


where A is the cross-sectional area of the now, C and k are empirieally-derived parameters, and u 
is the local flow velocity. By conservation, 
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Solving for C, the flow resistance coefficient, yields: 
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(3) 


Given measurements of cross-sectional area, thickness, and slope, as well as a value for k, this 
equation can be used to calculate relative changes in rheology from some initial point. The typica 
value k =0.5, employed in numerous previous studies, is adopted as a mathematical convenience. 
However, initial attempts to estimate k imply that the value for the Generals slide is somewhere m 
this vicinity. C„/Co values are given in Table 1 under “Modeling results”. Initial results from the 
application of this model are discussed in Bulmer et al. (Appendix 4). 
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Table la. Generals Slide 


Minimum Case 


Modeling Data 


MADISON DATA (GENERALS SLIDE) AS OF 3/27/00 MINIMUM 

Transect downstream position (m) elevation (m) trended slope (deg) width (m) weighted av. thickness (min) (m) 


AA 

473 

438 

BB 

535 

421 

CC 

627 

402 

Mad5 

808 

374 

Mad4 

905 

352 

Mad3 

981 

344 

Mad1 

1084 

323 

Mad2 

1173 

311 

E 

1361 

294 

DD 

1404 

290 


15.17 32.51 2.72 

14.05 28.01 1-83 

12.51 31.12 2.85 

9.9 44.68 1.62 

8.75 21.54 2.77 

7.98 33.69 3.12 

7.1 25.09 2.65 

6.5 28.93 3.39 

5.75 17.17 2.39 

5.67 24.55 2.23 


Used polynomial 3-order fit to calculate slope 
y=-5.6717E-8x3 + 0.0002544x2 - 0.472968x + 613.03 
R2 = 99.8107 

Based on horizontal flow surface; thicknesses are from weighted average depths (Area/width) 


Modeling 

Area results 

av. A (min) (m2) Cn/Co 


88.5 
51.3 

88.6 

72.5 

59.6 
105.2 

66.6 

98.1 

41.1 
54.7 


1 

0.475209 

0.867806 

0.354173 

0.485827 

0.663468 

0.399199 

0.567643 

0.229254 

0.243351 



Table 1b. Generals Slide 


Maximum Case 


Modeling Data 


MADISON DATA (GENERALS SLIDE) AS OF 3/27/00 

Transect downstream position (m) elevation (m) trended slope (deg) 



473 

BB 

535 

CC 

627 

Mad5 

808 

Mad4 

905 

Mad3 

981 

Madi 

1084 

Mad2 

1173 

E 

1361 

DD 

1404 


438 

15.17 

421 

14.05 

402 

12.51 

374 

9.9 

352 

8.75 

344 

7.98 

323 

7.1 

311 

6.5 

294 

5.75 

290 

5.67 


M/\XIMUM 


Area 


dth (max) (i weighted av. thickness (max) (m) av. A (max) (m2) 


49.02 

41.41 

53.61 

66.1 

21.54 

37 

34.72 

39.84 

18.69 

25.28 


3.81 
3.07 
3.05 
1.91 
2.77 
4.54 

3 

4.15 

2.82 
3.04 


186.6 

127.1 

163.3 

126.4 
59.6 
167.9 

104.4 

165.2 
52.8 

77 


Used polynomial 3-order fit to calculate slope 
y=-5.6717E-8x3 + 0.0002544x2 - 0.472968X + 613.03 
R2 = 99.8107 


Modeling 

results 

Cn/Co 

1 

0.616944 

0.619889 

0.271081 

0.238859 

0.599651 

0.278187 

0.443358 

0.150739 

0.193514 


k= 0.5 


Based on tilted flow surface; thicknesses are 


from weighted average depths (Area/width) 



TABLE 2. GENERALS DEBRIS 


FLOW: TRANSECT NOTES TABLE 


Transect 

AA 

BB 

CC 

Mad5 

Mad4 

Mad3 

Madl 

Mad2 

E 

DD 


Notes 

Just belowjmiction of two forks. Scoured bed. High “splash" on nver left. 

Scoured bedrock bed. 

Scoured bedrock bed. 

Flow separation, possible mass loss above this point. 

, , , a ^1 >jn clear field evidence defining right flow edge. 

Before/atop boulder dam. Well^constrained bedrock channel. No clear f.eld e 

Well-constrained bedrock channel. Definition of right flow edge problematic. 

tree-; Verv good flow depth/edge evidence. 

Well-constrained bedrock channel. Channel turns, emerges from trees. V ry g 

Deep incised V shaped clay-rich soil channel. Very good flow depth/edge evidence. 

V shaped clay-rich soil channel. Some evidence for deposition. Very good flow depth/edge evtdence. 

V shaped clay-rich soil channel. Some evidence for deposition. 
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TOPOGRAPHIC DATA FOR DEBRIS FLOWS - IMPLICATIONS FOR PLANETARY MODELING 
STUDIES. ’Bulmer, M.H., ^Peitersen’M.N., ^Bamouin-Jha,O.S., 'Johnston, A.K., and M. Bourke, Center^for Earth 
and Planetary Studies, Smithsonian Institution, Washington DC 20560-0315 (mbulmer@ceps..nasm.edu), Proxemy 
Research. 20528 Farcroft Lane, Laytonsville, MD 20882, 'The Johns Hopkins University Applied Physics Labora- 
tory, Laurel, MD 20723-6099. 

Introduction: Debris flows and other forms of mass feldspathic rocks [3]. Surface weathering has turned this 
movement are significant geomorphic agents on the ter- rock into loose friable soils with substantial qu^z s^d 
restrial planets. Viking images revealed deposits that and clay sitting on top of the coherent bedrock. Collu- 
have been interpreted to be similar to terrestrial debris vium formed from soil creep and slope wash overlays 
flows [1]. The term debris flow is used to denote a vis- some of the steeper slopes (slope >10°). Meters of 
cous to highly fluid form of rapid mass movement of weathered residual materials overlie bedrock at the base 
granular solids, (vegetation and air) with flow properties of steep slopes. .... 

that vary with water content, sediment size and sorting of To investigate the emplacement mechanism of debris 
particles. We have obtained high-precision topographic flows in the region we collected GPS topographic 
profiles and performed a sedimentological study over a profiles, and rock/soil samples on the 'Generals slide 
debris flow near Graves Mill (38.24°N, 78.23°W) in [Fig. 1]. During the topographic survey, slope angles, 
Madison County, Virginia. Use of an empirical model in flood-lines and boulder positions were also measured, 
conjunction with these data constrain to first order the The geometry of the channel both along and across the 
flow behavior during the emplacement of the Graves Mill stream profile were characterized using one long transect 
debris flows, and provide a method to interpret debris from the debris fan up to the source region, and ten 

across channel profiles where ever it was possible to 
match up high-water marks on the left and right channel 
banks. 

Field results: Geomorphic evidence indicates that more 
than one pulse formed the Generals slide debris flow 
deposits. The actual number is unknown, but there are 
several drainage channels which feed into the main 
channel [Fig 1.]. The distance from the top of channel 'a' 
(dl) of the 'Generals slide' to the edge of the treeline (d3) 
is 919 m and from the top of channel 'b' (d2) to the same 
location (d3) is 925 m. The elevation change over this 
length is 329.3 m. There are many breaks of slope and 
good evidence for superelevation on the outside of 
channel bends. At the source area to channel 'a' the 
failure scarp is only 14.7 m wide and 0.7 m to 1 .0 m deep 
on a 48° slope. This initial failure appears to have been a 
shallow translational soil and rock slide. The channel the 
debris flow traveled down contains a range of coarse- 
sized materials with boulders up to 10 m on a side, some 
of which are from a prior event but are mantled by clasts 
from the 1995 event. Once the debris flow emerged from 
the channel near the treeline and out onto the fields [Fig. 
1] it deposited a range of material sizes from blocks and 
boulders down to clay and silt. Based on measurements 
on an air photo the main deposit L (d3 to d4) is 3 1 1 m 
long and 93 m at its widest W [Fig. 1]. It covered an area 
0.02 km'. Eyewitness accounts estimate the thickness T 
of the deposits to have been approximately 1 m. Using L, 
W and T we have derived a first-order volume estimate 
of 21,960 m' of deposited material. This is significantly 
greater than that which could have been derived from the 
source areas 'a' and 'b' alone. 


flow emplacement on Mars [e.g. 1]. 



Field Site: An intense storm on June 27, 1995 centered 
over an area of about 7 5 km' triggered hundreds of rock, 
debris and soils slides from the steep hillsides of the She- 
nandoah Mountains [2]. Most of these slope failures 
transformed into debris flows [Fig. 1]. During approxi- 
mately sixteen hours as much as 770 mm of rain fell in 
the area of maximum storm intensity and probably -640 
mm fell within a five-hour period over small areas [2]. 
The area is underlain by mostly granitic quartzo- 


A v>T> rL . A 
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topographic data for a debris FLOW: Bulmer, M.H., al. 
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The Empirical Model; Simple empirical volumetnc 
now rate models provide preliminary insight into the 
behavior and rheology of the Generals slide. Such models 
have been used to study diverse mass movements of 
geologic materials, including lahars [4], mud flows [5], 
glaciers [6], and lava flows [e.g. 7,8]. In this study, we 
use a very general form for flow rate 

gh sin0^* 


that will further constrain the flow dynamics as expressed 
in the empirical model used here. 


Q = uA = 


(1) 


Q _ Fr„ _ 


^^0 


where A is the cross-sectional area of the flow, g is grav- 
ity 0 is the underlying slope, h is the flow thickness, and 
C and k are empirically derived parameters. The vanable 
C = 1/Fr when k = 0.5. The Froude number Fr relates 
inertial to gravitational forces and thereby C is indicative 
of all the other forces relative to gravity that influence the 

downhill progress of the debris flow. ^ ^ 

We can use the morphometry of die debris flows t 
estimate C. As in other studies of this type [e.g. we 
maintain k = 0.5 to ensure a dimensionless C. In addi- 
tion, we assume a constant Q. Using our morphometric 
data set, we estimate changes in C at any point along 
the debris channel relative to a reference point xo'. 

2 

hg sin9o ^2) 

A„J^„sin0„ 

Inferences about processs along the channel can be made 
by comparing C to morphologic features downstream, as 
well as to parameters found for other debris flows. 

Results of model application/discussion: Application of 
equation (2) revealed large changes in Fr, as shown in 
figure 2. Initially decreasing Fr are consistent with a 
debris flow becoming less fluid possibly due to entr^ 
ment of eroded material. Large variation at around 1000 
m occurs near rapid slope and cross-sectional area 
changes. High values of Fr from the 1000 m point to the 
flow’s end are most likely overestimates in C caused by 
the constant Q assumption. A drop in Q due to loss of 
-volume in the depositional zone is expected, and wou 
result in increasing Fr as shown by figure 2. The model 
estimates for flow speed along the channel will be tested 
against superelevation data as well as observed blocks 
sfzes on the flow margins. Such comparisons will better 
constrain the C, Fr and Q along the Generals s ide 
thereby improving our preliminary understanding of the 
flow rheology as it moves downslope. 

Planetary implications: These preliminary analyses 
indicate how high resolution topographic data may pro- 
vide a first step towards constraining the dynamics and 
rheology of a debris flows. The block size data, the su- 
perelevation analysis and the sedimentary analyses along 
the channel of the flow provide additional information 




0 250 500 750 1000 1250 1500 

Distance downstream (m) 

Fig. 2. Variations in elevation, cross-sectional area A and 
Fr^/Fr^ with distance downstream. 

In addition, this information can be linked to more 
sophisticated models that allow estimating the amounts of 
fluid and solid [e.g. 9] in a debris flow, information of 
particular importance to understanding 
evolution of water reservoirs on Mars. On M^s, 
topographic data can be obtained from MOLA while 
sedimentary data may be obtained through thermal inertia 
obtained at a debris flow using TES data, as well as hig 
resolution MOC imagery. 

References: [1] Carr, M.H.. (1996) Water on 
Wieczorek. G.F., et al.. (1996) USGS Open Fi e Rep^ 96- T 
[3] Allen, R.M. (1963) Virginia Div. Mm. 4 

Weir G.J., (1982) New Zealand J. Set. 25. 

Gol’din, B.’and L.S. Lyubashevskiy (1966) Trudy l^kr. NIC 
60, 73-75. [6] Paterson. W.S.B. (1969) Physics of'Glac.erj 
250 pp. [71 Baloga, S.M, et al., (1995), J. Ceophys. Res. m, 
24 509-24519. [8] Baloga. S.M., et al. (1998) J. Geophys^ Res_ 
idl 5133-51421. [9] Iverson, R. M. (1997) Rev. Geophs. 35, 
3. 245-297. 
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Investigating the physical factors controlling a debris flow - Implications for landscape 
evolution. Bamouin-Jha, O., M. Bulmer, M. Peitersen, M. Bourke and A. Johnston. 

Detailed field work of a debris flow near Graves Mill (38.24N, 78.23W) in Madison County 
Virginia is combined with theoretical modeling to better constrain the moa important physical 
factors controlling debris flow mechanics. Such work will allow improved characterration of 
Lbris flows and Lir effect on landscape evolution, and help interpret debris flow emplacement 
on other terrestrial planets (e.g., Mars). 

T vwia O'? 100 ^ triocTprfsfl 3 . siihstantial number of debris flows off the steep 

The now (hemafter denoted as the 

“Generals Slide”) extends from an obvious source area (a translational sli^ sc^ possessing 

4 " slope) to a depositional fan approximately 1230 m “ 

m hiah resolution GPS topography both along and across the channel, (2) extens ve 

7" 7:tix 

data at various bends in the channel; (4) sedimentary profile^ (5) s^d|inentary 
and boulders comprising the debris flow; (6) cohesion and internal friction of the 
debris flow’s source, and (7) the plastic and liquid limits of the matrix. 

In a previous study the topographic data was input into a simple empirical volumetric flow rate 
Idel preliminary insight into the behavior of the General slide. The resuUs 

indicate qualitative coLlations between the flow dynamics as expressed by this model with run- 
off input local changes in slope, and sediment loss and gain. In the current work, ngorous^^ 
test these model results against quantitative information obtained from our fieldwork data for the 
debits florspe^^^^^ theology, and sedimentary profiles. We also assess more sophisticated 
ILs flow models that include both the measured topography and the rheological characteristics 

of the flow. 
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Topographic Data for Debris Flows - Implications for Planetary Modeling Judies. 

'Bolmer, M.H., =Pcitersen, M.N., =Bamouin-Jha.O.S., 'Johnson, A.K., and M. ‘Bourke, 
'center for Earth and Planetary Studies, Smithsonian Institution, Washington DC 20560- 
0315 (n.buln,er@ceps.nasm,edu), Wemy Reseamh, 20528 Farcroft Lane, Laytonsville, 
MD 20882, ’The Johns Hopkins University Applied Physics Laboratory, Laurel, MD 

20723-6099. 


Abstract. Debris flows and ote forms of mass movement are significant geomorphic 

interpreted to be similar to terrestrial debris flows [Cnrr, 1996], However, many 
questions about the origin of these deposits can only be addmssed through high-resolution 
imaging, topography and sedimentologiea. data. The term debris flow is used to denote a 
viscous to highly fluid form of rapid mass movement of granular solids, (vegetation and 
air) with flow properties that vary with water content, sediment size and sorting 
particles We have obtained high-preeision topographic profiles and performed a 

lento, ogica. study over a debris flow near craves Ml, 1^ 

Madison County, Virgima. Use of an empirical model in coniunetion with these data 
constrain to first order the flow behavior during the emplacement of the Graves Mill 
debris flows, and provide a method to interpret debris flow emplacement on Mars [Cure, 

1996]. 
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Introduction 

2 

An intense storm on June 27, 1995 centeted over an area of about 75 fan 
triggered hundreds of rock, debris and soils slides from the steep hillsides of the 
Shenandoah Mountains [If/cczoret el ai. 1996], Most of these slope failures transformed 
into debris flows (Figure 1), During approximately sixteen hours as much as 770 mm of 
rain fell in the area of maximum storm intensity and probably -640 mm fell wtthtn a five- 
hour period over small areas [Wieczorek el ai. 1996], Numerous geologic studies have 
been carried out on debris flows [eg., Cosla and Wieczorek. 1987; Iverson. 1997], but 
many questions remain that need to be addressed using detailed field studies before 
possible debris flow features on Mars can be adequately examined. Historical debris 
flows have been reported all along the Appalachian Mountains, and stratigraphic 


evidence indicates that debris flows have been important geomorphic slope processes 
through the Quaternary [Mills. 1982; Kochel and Johnston. 1984], An inventory of debris 
flow and flooding effects from the 1995 storm in Madison county was produced by the 
uses [Wieczorek et ai. 1995, 1996; Morgan el ai. 1997], The area is underlain by 
mostly granitic quartzo-feldspathic rocks [Allen. 1963], Surface weathering has turned 
this rock into loose friable soils with substantial quartz sand and clay sitting on top of the 
coherent bedrock. Colluvium formed from soil creep and slope wash overlays some of the 
steeper slopes (slope >10"). Meters of weathered residual materials overlie bedrock at the 

base of steep slopes. 

We chose to investigate the emplacement mechanism of one debris flow from the 
1995 storm known as the “Generals slide" (Figure 1). We collected topographic profiles 
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using global positioning (GPS) equipment, and rock/soil samples from both the “Generals 
slide" deposits and along the stream channel. During the topographic survey, slope 

angles, flood-lines and boulder positions were also measured. 

In this paper, we consider the use of an empirical model in conjunction with the 
geomorphic observations and stratigraphic studies, to constrain to firs, order ttte 
emplacement of the “Generals slide” debris flow and hence a method to interpret debris 
flow emplacement on Mars. We firs, describe flte techniques used to obtain the 
topographic dataset from the source area all the way out to the depositional fan of the 
“Generals slide”. Profiles both along the channel and fan as well as at selected locati 
across the channel. The results from the fieldwork surveys are then described. These data 
are then used in the flow rate model. We val.date flte model-derived outputs by 
comparing them with our geomorphological and sedimentological observations. Finally, 
we assess flre impact of these results to studies of possible Martian debris flows. 


Surface Topography Data Collection 

The main channel of the “General slide” was surveyed using a Trimble Total 
Station 4800, which consists of two GPS receivers that work in hmdem. This equipment 
uses canier-phase differential processing for 1 cm horizonuil real-time accuracy. One 

tripod at a known location during the survey (Figure 2). The 

of 1-2 m using a Trimble ProXRS 


receiver remains fixed on a 

location of this receiver was determined to an accuracy 

lite-based differential correction. This base station was then used as 


GPS unit using satellite- 
the reference point for all subsequent surveys. The roving 


Trimble 4800 receiver used for 
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data collection is mounted on a pole 1.8 m in 
determine positions even in the small spaces 


length. This pole allows the survey team to 
between blocks, so the topography can be 


adequately eharacterized. A radio link, provided by a 25 W UHF radio at the base station, 
is constantly maintained for real-time differential correction. While conducting transects 
up the main channel, a 3-D point was collected approximately every 3 m. This spacing 
was determined in the field to adequately represent the topographic character of the 
channel profile and met the resolution requirements of the empirical model. For the 
transeet up the channel (transect 1) we obtained 1 cm horizontal and 2 cm vertical 
accuracy per point to about 1000 m up from the lower end of the channel, nntil we 
encountered difficulties in satellite coverage due to screening by the tree canopy. Beyond 
that point, the accuracy was only about 20-40 cm. We successfully collected data up to 

the top of the main channel, about 1600m. 

In addition to the main channel transect we also measured 10 cross-channel 

profiles (Figure 2) using a laser rangefinder (Laser Atlanta Advantage CL) connected to 
Uie Trimble ProXRS. A 3-D point was collected every 3-6 m with an accuracy of 10-15 
cm. The position of each cross-channel profile relative to the main channel was 
determined to an accuracy of about 1 m by locating points that had been surveyed 
pmviously using the Trimble 4800. Using available air photos and field observations, the 
attempt was made to conduct cross-channel transects where a measurement of depth for 
the 1995 event could be determined as well as an across-channel profile. Geomorphic 
signatures of the 1995 event that we used frequently in determining debris flow depths 
were locations where tree, soil, and rock materials had been pushed up against the base of 
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,r.e trunks on the edge of the channel. Measurements of minimum depths were also 
colleeted where fines from the 1995 even, had been deposited on top of large boulders. At 
some locations (transects BB. CC, Mad4) it was difficult to determine the edge of the 
1995 event. At these locations we extended the profile out to the first 
geomorphic indicator for the absence of the ground having been disturbed by the 1995 
event. Using this procedure we attempted to obmin the maximum possible width 

measurements. 


Field results 

Geomorphic evidence indicates that more than one pulse formed the Generals 
slide debris flow deposits and that these pulses were likely ttansitional between 
hyperconcentrated flow and debris flow. The actual number of pulses is unknown, but 
there are several drainage channels, which feed into the main channel (Figure 1). The 
disumce from the top of channel 'a' (d,) of the ■Generals slide’ to the edge of the treeline 
(ds) is 919 m and from the top of channel 'b (dt) to the same location (d,) is 925 m. The 
elevation change over this length is 329.3 m. There are many breaks of slope and good 
evidence for superelevation on the outside of channel bends. At the source area to channel 
■a' the failure scarp is only 14.7 m wide and 0.7 m to 1.0 m deep on a 48" slope. This 
initial failure appears to have been a shallow translational soil and rock slide. The channel 
the debris flow traveled down contains a mnge of coarse-sized materials with boulders up 
to 10 m on a side, some of which are from a prior event but are mantled by clasts from the 
1995 event. Once the debris flow emerged from the channel near the treeline and out onto 
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the fields (Figure 1) it deposited a range of material sizes from blocks and boulders down 
to clay and silt. Based on measurements on an air photo the main deposit L (d 3 to d 4 ) is 
311m long and 93 m at its widest W (Figure 1). It covered an area 0.02 kml Eyewitness 
accounts estimate the thickness T of the deposits to have been approximately 1 m. Using 
L, W and T we have derived a first-order volume estimate of 21,960 m^ of deposited 
material. This is significantly greater than that which could have been derived from the 
source areas 'a' and 'b' alone. 


The Empirical Model 

Simple empirical volumetric flow rate models provide preliminary insight into the 
behavior and rheology of the Generals slide. Such models have been used to study diverse 
mass movements of geologic materials, including lahars [JVier, 1982], mudflows [Go/ ’din 
and Lyubashevskiy. 1966], glaciers [Paterson, 1969], and lava flows [e.g. Baloga et al., 
1995, 1998]. In this study, we use a very general form of the flow rate (Q); 

Q = uA = [gh sin 6/ C]'‘ A (^) 

where A is the cross-sectional area of the flow, u is the velocity, g is gravity, ^is the 
underlying slope, h is the flow thickness, and C and k are empirically derived parameters. 
Assuming a steady state and a constant flux. 

Qo = [gho sin 9o/ Co] An = [gh„ sin 9„ / C„] * = Qn (2) 
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Solving for C, the flow resistance coefficient, yields 


[C„/ Co] = [AJ Ao] / [(hn sin 0„ )/ (ho sin ^o)] 

Given the morphometry of the debris flow, and a known value of k, this equation can be 
used to calculate relative changes in rheology (as expressed by variations m C) from some 
arbitrary initial reference point.. The typical value k = 0.5, employed in numerous 
previous studies type [e.g. Bologa et ai, 1998], is adopted as a mathematical 
convenience. However, initial attempts to estimate k imply that the value for the 
Generals slide is somewhere in this vicinity. 

Results of model application/discussion 

Application of equation (3) revealed moderate changes in C, as shown in Figure 3. A 
decrease in C with distance implies increased fluidity, consistent with the debulking of 
the debris flow by additional surface runoff. Small-scale topographic effects may cause 
deviations from this trend. A comparison of the flow profile to the modeling results 
(Figure 3) reveals positional correlations between repeated local slope changes and 
variations in C, in some cases. Surface roughness, nature of the bed, and flow/bed 
interactions (erosion, deposition, etc.) as well as local variations in volume (loss from 
overbanking, increase from tributary rills), may be responsible for the remaining 
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variations. These factors, which, can be best constrained by field observations, will be 
subjects of future study. 

Planetary implications 

These preliminary analyses indicate how high resolution topographic data may 
provide a first step towards constraining the dynamics and rheology of a debris flows. 
Sufficiently accurate elevations, traditionally lacking for planetary features, are currently 
available for Mars via the MOLA data. However, topographic data alone is insufficient 
to adequately understand the dynamics of the flow. The block size data, the 
superelevation analysis and the sedimentary analyses along the channel of the flow, as 
well as detailed field observations, provide the additional information that will further 
constrain the flow dynamics as expressed in the empirical model used here. Such 
information will largely be lacking for planetary debris flows; however, judicious use of 
additional data sets such as TES and MOC, as well as comparison to appropriate 
terrestrial analogs, may be used to constrain Martian flows. 

Our field data shows the complex processes involved in the debris flows. In 
future work, this information will be linked to more sophisticated models that allow 
estimating the amounts of fluid and solid [e.g. Iverson, 1997] in a debris flow, 
information of particular importance to understanding the size and evolution of water 
reservoirs on Mars. This work provides a foundation for analyzing existing MOC, 
MOLA, TES and Viking data over proposed Martian debris flows. 
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FIGURE CAPTIONS: 


Figure 1. Subscene from air photo of the Graves Mill area, taken 22 August 1995, several 
months after the debris flow event. Depositional fan of Generals slide is still clearly 
evident. 

Figure 2. Field survey transects, overlaid on air photo from 24 March 1997. 

Figure 3. Results of application of empirical model to Generals slide data. Relative 
change in C is overlaid on a plot of transect 1 . 
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